Two Scenarios of the Quantum Critical Point 
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Two different scenarios of the quantum critical point (QCP), a zero-temperature instability of 
the Landau state, related to the divergence of the effective mass, are investigated. Flaws of the 
standard scenario of the QCP, where this divergence is attributed to the occurrence of some second- 
order phase transition, are demonstrated. Salient features of a different topological scenario of the 
QCP, associated with the emergence of bifurcation points in equation e(p) — fj, that ordinarily 
determines the Fermi momentum, are analyzed. The topological scenario of the QCP is applied to 
three-dimensional (3D) Fermi liquids with an attractive current-current interaction. 
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A statement that the Landau quasiparticle picture 
breaks down at points of second-order phase transitions 
has become a truism. The violation of this picture is at- 
tributed to vanishing of the quasiparticle weight z in the 
single-particle state since the analysis of a long wave- 
length instability in the S = 1 particle-hole channel, 
performed more than forty years ago by Doniach and 
Engelsbergi and refined later by Dyugaev.^ In nonsuper- 
fluid Fermi systems, the z-factor is determined by the 
formula z = [1 — (9E(p, £)/9e)Q]~^ where the subscript 
indicates that the respective derivative of the mass op- 
erator E is evaluated at the Fermi surface. This factor 
enters a textbook formula 
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for the ratio M*/M of the effective mass M* to the 
mass M of a free particle. As seen from this formula, 
where = p^/2M, the effective mass diverges at a 
critical density pc, where z vanishes provided the sum 
1 -t- (91] (p, e = 0)/9ep)p has a positive and finite value at 
this point. Nowadays, when studying critical fluctuations 
of arbitrary wave-lengths k < 2pp has become popular, 
this restriction is often assumed to be met without stipu- 
lations. E.g. a standard scenario of the quantum critical 
point (QCP) where M* diverges is formulated as follows: 
in the vicinity of an impending second-order phase tran- 
sition, "quasiparticles get heavy and die".-'- 

However, as seen from Eq.(IT]), M* may diverge not 
only at the points of the second-order phase transi- 
tions, but also at a critical density poo, where the sum 

1 + changes its sign. Furthermore, we will 

demonstrate that except for the case of the ferromagnetic 
instability,! M* cannot diverge at pc without violation of 
stability conditions. 

In what follows we restrict ourselves to one-component 
three-dimensional (3D) homogeneous Fermi liquids 
where the particle momentum is conserved, and the Lan- 
dau equation, connecting the quasiparticles group veloc- 
ity de/dp to their momentum distribution n(p) in terms 



of the interaction function /, has the form^i^ii 



dp 



(2) 



Setting here p = pp and introducing the notation vp 
{de{p)/dp)Q —Pp/M*, one obtains 



PF A PfM 



implying that 
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Hereafter we employ notations of Fermi liquid (FL) the- 
ory where /i is the first harmonic of the interaction func- 
tion f{e) = z'^T'^{pp,pp;9), with T'^ being the w-limit 
of the scattering amplitude F of two particles, whose en- 
ergies and incoming momenta Pi,P2 he on the Fermi 
surface, with cos6 = pi • P2/p'p, while the 4-momentum 
transfer (q, w) approaches zero, such that q/uj — s- 0. 

It is instructive to rewrite Eq.(|3]) in terms of the k- 
limit of the dimensionless scattering amplitude i/F'^ — 
A + B(Ti(T2 where ly — z'^ppM* /tt'^ is the quasiparticle 
density of states. Simple algebra then yields 



M 

AI* 



1 



(5) 



This formula stems from Eq.ig]) and relation^-^" Ai = 
$i/(l + ^), where $ is the spin-independent part of the 
product i^F" . Thus at the density poo where the effective 
mass diverges one has 

^l(/5oo)=3, $i(poo)=oo. (6) 

In the following we focus on critical density fluctua- 
tions with fcc 7^ 0, addressed in Ref. S. First we no- 
tice that there is a strong dependence of the amplitude 
rQ/3,75(pi, P2, k, a; = 0) on the momentum transfer k 
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FIG. 1: Singular contributions to the scattering amplitude in 
the vicinity of the second-order phase transition. 



close to the critical momentum kc, specifying the spec- 
trum of density fluctuations, that stems from the asym- 
metry of r with respect to the interchange of momenta 
and spins of colliding particles.^ In this case, upon ne- 
glecting regular components one finds^, (see Fig.[T]), 

A(pi,p2,k,w 0;p^ Pc) = -L»(k) + ii:>(pi-p2 + k) , 



with 



D{k ^kc,uj = 0) 



+ {k- kj^ 



(7) 
(8) 



the correlation length ^(p) diverging at p — pc- 

Within the quasiboson approximation®, the derivative 
{d'S{p,e)/de)Q diverges at p ^ Pc as ^(p), while the 
derivative e)/9ep)p remains finite at any kc- If 

these results were correct, then the densities pac and pc 
would coincide, in agreement with the standard scenario 
of the QCP. However, calculations of harmonics Ai;{p) of 
the amphtude A{pp,pF, cos9) from Eqs.l?]) and ^ yield 



A I \ '^^ci{p) , , . 



g— 2 — cos6io. 

2 Pf 



(9) 

We see that the sign of Ai{p pc), coinciding with 
that of cos^o = 1 ~ ^c/^Pfj turns out to be negative 
at kc > pfV^- According to Eq.fl]) this implies that at 
the point of the second-order phase transition, the ratio 
M*{pc)/M < 1. Thus we infer that at kc > pfV^, the 
densities pc and p^c cannot coincide. In its turn, this 
implies that vanishing of the z-factor at pc is compen- 
sated by the divergence of the derivative {d'E{p,e)/dep)o 
at this point, otherwise Eq.([T|) fails. 

To verify this assertion let us write down a funda- 
mental FL relation^ between the k- and a;-limits of 
the vertex T that has the symbolic form T'^ = T'^ -t- 
(r'=((G2)'= - (G'2)'^)r") where external brackets mean 
integration and summation over all intermediate mo- 
memtum and spin variables. In dealing with the bare 
vertex T° — p the extended form of this relation is 

dG-\p) _ dG-\p) p 



dp 



de M 
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dG-\q) q d^q 
dqa M (2^ 
(10) 



In writing this equation the Pitaevskii identities^ 

r^-(p) = -^^^, r-(p) = ^^^P (11) 

are employed. Upon inserting the FL formula 

(G^q))' ^{G\q))^ -2n^t^S{e)6{p^PF) (12) 

Pf 

into Eq. (fTO|) and the standard replacement of the spin- 
independent part of i'T'' by A, after some algebra we are 
led to equation 
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Remembering that —dG ^(p, e)/9p — z ^de{p)/dp one 
arrives^ at Eq.®. On the other hand, as seen from 
Eg. pT])) . at kc > PfV^ where Ai < 0, the derivative 
(d'E{p,e)/dep)^ does diverge at the same density, as the 
derivative (9S(p, e)/9e)g, in contrast to the result.— To 
correct the defect of the quasiboson approximation®, the 
spin-independent part of the scattering amplitude F en- 
tering the formulas® for the derivatives of the mass op- 
erator E should be replaced by that of the amplitude F" 
(for details, see Ref. y)- 

We will immediately see that at finite kc < pfV^, 
vanishing of the z-factor is incompatible with the di- 
vergence of M* as well. Indeed, as seen from Eq.®, 
the harmonics Aq{pc) and Ai{pc) are related to each 
other by equation Aq{pc) = Ai{pc)/{3cos9o). If M*{pc) 
were infinite, then according to Eq.®, Ai{pc) would 
equal 3, and Ao{pc) = 1/ cos^o. However, the quantity 
Aq = $0/(1 + 'i'o) cannot be in excess of 1, otherwise the 
Pomeranchuck stability condition^ii $0 > — 1 is violated, 
and the compressibility turns out to be negative. Thus 
the QCP cannot be reached without the violation of the 
stability condition. If so, approaching the QCP, the sys- 
tem undergoes a first-order phase transition, as in the 
case of 3D liquid '^He. 

The finiteness of M* at the points of vanishing of 
the z-factor requires an alternative explanation, (see e.g. 
Ref. [^), of the logarithmic enhancement of the specific 
heat C(T), observed in many heavy fermion metal&iS, and 
attributed to contributions of critical fluctuations. 

Let us now turn to the analysis of another opportunity 
for the occurrence of the QCP, addressed first in micro- 
scopic calculations of the single-particle spectrum of 3D 
electron gas.^^'^^ It is associated with the change of the 

sign of the sum 1 -|- (^^-^§r^^ at poo ^ Pc- In this case, 

the z-factor keeps its finite value, and hence, the quasi- 
particle picture holds on both the sides of the QCP. In 
standard Landau theory, equation 



e(p, T = 0) = p 



(14) 



with being the chemical potential, has the single root, 
determining the Fermi momentum pF- Suppose, at a 
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FIG. 2: Single-particle spectrum e(p) of a homogeneous 2D 
electron gas in units of Sp — p%/2M, evaluated^** at T = 
for different values of the dimensionless parameter — 
\/2Me^ Ipf. 



critical coupling constant gr, a bifurcation in Eq. (jl4p 
emerges, then beyond the critical point, at <? > gr, this 
equation acquires, at least, two new roots that triggers 
a topological phase transition. ^'^ In many-body theory, 
equation, determining critical points of the topological 
phase transitions, has the form 



(15) 



Significantly, terms, proportional to elne, existing in the 
mass operator S of marginal Fermi liquids, do not enter 
this equation. 

The bifurcation pb in Egs. fH)) and can emerge 
at any point of momentum space. If pb coincides with 
the Fermi momentum pp, then at the critical density 

the sum 1 + (^^^^^^ vanishes, and one arrives at the 

topological quantum critical point. In connection with 
this scenario, it is instructive to trace the evolution of 
the group velocity vp — {de/dp)^ versus the first har- 
monic /i. As follows from Eq.([3]), vp keeps its positive 
sign as long as = ppAIfi/Sn"^ < 3, and the Lan- 
dau state with the quasiparticle momentum distribution 
np{p) = d{p — pp) remains stable. However, at > 3, 
the sign of vp changes, and the Landau state is neces- 
sarily rearranged. This conclusion is in agreement with 
results of microscopic calculations of the single-particle 
spectrum e{p,T ~ 0) of 2D electron gas,^'' shown in 
Fig. [21 As seen, the sign oi vp holds until the dimen- 
sionless parameter attains a critical value r^^ — 7.0. 
At greater r^, the derivative {de{p) / dp)^^, evaluated with 
the momentum distibution np{p), becomes negative, and 
the Landau state loses its stability, since the curve e{p) 
crosses the Fermi level more than one time. 

At T = 0, two types of the topological transitions are 
knownii^ One of them, giving rise to the multi-connected 
Fermi surface, was uncoveredi^ and studied later— li^iiiii^ 
in calculations of the single-particle spectrum e{p) on the 
base of Eq. with the interaction function /(fc), hav- 
ing no singularities at fc ^ 0. In this case, beyond the 
QCP, Eg. pi)) has three roots pi < P2 < P3, i-e. the 



curve e(p) crosses the Fermi level three times, and occu- 
pation numbers are: ■n{p) = 1 at p < pi, n{p) = at 
Pi < p < p2, while at p2 < p < P3, once again n{p) — 1, 
and aA, p > P3, n{p) — 0. As the coupling constant g in- 
creases, the number of the roots of Eq. (fT4|) rapidly grows, 
however, their number remains countable at any g > gx- 
The situation changes in Fermi liquids with a singular 
attractive long-range current-current term 



r°(pi,P2,k,Lu = o) 



P1P2 - (pik)(p2k)/fc^ 



(16) 

since in these systems, e.g. in dense quark-gluon plasma, 
solutions with the multi-connected Fermi surface are un- 
stable. Indeed, the group velocity de{p)/dp evaluated 
with np{p) = 9{p — pp) from Eq. 12]) has the form 
de{p)/dp ^ Pp/M - g\n{2pp/\pp - p\), implying that 
Eq.(fT4|) has three different roots pi,p2,?'3, corresponding 
to the Fermi surface, having three sheets at any 5 > 0.— 
However, at the next iteration step, the new Fermi sur- 
face has already five sheets, the Fermi surface has seven 
sheets and so onJ^ With increasing the number of it- 
erations, the distance between neighbour sheets rapidly 
shrinks. In this situation, a minute elevation of tem- 
perature renders the momentum distribution n{p,T) a 
smooth T-independent function n*(p), different from 
and 1 in a domain C between the sheets. In this case, the 
ground-state stability condition. 



6E = J {e{p) 



fi) 5n{p) 



d^p 
(2^ 
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(17) 



requiring the nonnegativity of the variation 5E of the 
ground state energy E at any admissable variation of 
n^{p), is met provided 



e(p) = 11- , peC . 



(18) 



As a result, we arrive at another type of the 
topological transitions, the so-called fermion 
condensatioi*i^i^'2iiS^i2^iHi2^, where the roots of 
Eq. (|14p form an uncountable set, called the fermion 
condensate (EC). Since the quasiparticle energy e{p) is 
nothing but the derivative of the ground state energy E 
with respect to the quasiparticle momentum distribution 
n{p), Eq. p^ can be rewritten as variational condition — 



6E 



At, pGC . 



5n{p) 

The EC Green function has the form 
1 - (p) [p) 



G{p,e) 



iS 



i6 



peC 



(19) 



(20) 



As seen, only the imaginary part of the EC Green func- 
tion differs from that of the ordinary EL Green function. 
This difference exhibits itself in a topological charge, given 
by the integral^^'^^ 



N= / Gip,s)diGip,e) 



dl 
27ri 



(21) 
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where integration is performed over a contour in com- 
plex energy plane, embracing the Fermi surface. For 
conventional Fermi liquids and systems with the multi- 
connected Fermi surface, the topological charge N is in- 
teger, while for the states with a FC, its value is half- 

For illustration of the phenomenon of fermion conden- 
sation, let us address dense quark-gluon plasma, on the 
Lifshitz phase diagram of which, as we have seen, there is 
no room for the conventional FL phasci^ Upon inserting 
into Eci.([2]) only leading divergent terms in the interac- 
tion function /, constructed from Eq. ()16|l . one finds 

0^1- A /in , ^ x,xieC,{22) 
J \x-xi\ dxi 

where dimensionless variables x = {pp — p)/2pp and A 
are introduced. A numerical solution of this equation will 
be found elsewhere. Here we simplify Eq. (j22p . replacing 
the kernel ln(l/|a; — xij) by ln(l/a;) provided a; > xi, and 
by ln(l/xi), otherwise, to obtain 

= 1 + An* (.t) In .T + A / Inxi — *^ dxi x,xi^C. 

J OXi 

X 

(23) 

An approximate solution of this equation is n,(a;) = 
x/Xm where Xm, — is determined from condition 

n{xm) — 1- We see that the range of the interval [0,Xm] 
of fermion condensation, adjacent to the Fermi surface, 
is exponentially small. 

Nontrivial smooth solutions n* (p) of Eq. (fl9|l exist even 
in weakly correlated Fermi systems. However, in these 
systems, the Pauli restriction n^(p) < 1 is violated, ren- 
dering such solutions meaningless. Even at the QCP, 
where the nonsingular interaction function / is already 
sufficiently strong, no consistent FC solutions n^{p) exist, 
satisfying the requirement n{p) < 1 wherever. These so- 
lutions emerge at a critical constant gpcj and at g > gpc, 
they win the contest with any other solutions. Thus on 
the Lifshitz phase diagram of systems with nonsingular 
repulsive interaction functions /, the standard FL phase 
occupies the interval g < gx, the phase with the multi- 
connected Fermi surface, the interval gr < g < gpc, 
while the phase with the FC exists at g > gpc- 

In dealing with the full (T — g) phase diagram of such 
systems we notice that the temperature evolution of the 
quasiparticle momentum distribution, associated at T = 
with the multi-connected Fermi surface, depends on the 
departure of the difference |e(p) — /i| from in the domain 
C. Its maximimu value £m determines a new energy scale 
em — (P /M*(0), where d is the average distance between 
the sheets of the Fermi surface that rapidly falls with 
the increase of the sheets mimheiJ^^ If temperature T 
attains values, comparable with Cm, then, as seen from 
the Landau formula n{p) = [1 + exp (e(p) — /z)/T)]~^, 
the distribution n{p) becomes a smooth function of p. 
Employing the FC notation n^{p) for this function, one 



finds that at T > e™ the spectrum 

e(p,T) = rini^^, pec (24) 

does coincide with the FC spectrum22. We infer that 
at T ~ em, a crossover from a state with the multi- 
connected Fermi surface to a state with the FC occurs. 
As a result, FL thermodynamics of the systems with 
the multi-connected Fermi surface completely alters at 
T ~ e™, since properties of systems with the FC resemble 
those of a gas of localized spins. Such a transition was 
recently observed in the heavy-fermion metal YbIr2Si2, 
transition temperature being merely 1 Ki^ 

Let us now turn to systems of fermions, interacting 
with a "foreign" bosonic mode, e.g. phonons or pho- 
tons. In the Frolich model,— aimed for the elucidation 
of electron and phonon spectra in solids, electrons share 
momentum with the lattice due to the electron-phonon 
interaction. The non-conservation of the electron mo- 
mentum results in the violation of the second of the re- 
lations PT|) . and Eq.(|3]) acquires the form 

M^^(nr^/ 1 A 
M* Pp \ S J 

where n = p/p. The departure of the ratio M/M* 
from the Landau value is well pronounced in the 
limit Cs << vp. For illustration, let us consider the 
weak coupling limit of the Frolich model, where the 
first harmonic Ai, evaluated from the phonon propaga- 
tor D{\Tpi — P2|,t^ = 0) equals due to its isotropy. If 
Eq.(|4]) were correct, then M* /M would equal 1. How- 
ever, this is not the case: M* /M = l+g'^ppM/2n^ where 
g is the electron-phonon coupling constant.-- Evidently, 
if the ratio vp/cg drops, then the departure from Eq. ^ 
falls due to weakening of the contribution of the pole of 
the boson propagator. Such a situation occurs just in the 
vicinity of the QCP, since at this point vp = pp/M* = 0. 

So far information on the QCP properties of Fermi 
liquids is extracted from measurements, carried out in 
2D electron gas, 2D liquid ■^He and heavy-fermion met- 
als. Here we restrict ourselves to several remarks, re- 
serving a more detailed analysis for a separate paper. 
Accurate measurements of the effective mass M* in di- 
lute 2D electron gas are made on (100)- and (lll)-silicon 
MOSFET's,^2i3ii32i33 In principle, the divergence of M*, 
observed in these experiments, can be associated with 
critical spin-density fluctuations. However, experimental 
data rules out a significant enhancement of the Stoner 
factor. In many theories, (see e.g. Ref. [H), the enhance- 
ment of M* is related to disorder effects. However, the 
effective masses, specifying the electron spectra of (100)- 
and (lll)-sihcon MOSFET's, where disorder is differ- 
ent, almost coincide with each other provided dimension- 
less parameters of the 2D Coulomb problem, are the 
samej2^ On the other hand, this coincidence that agrees 
with results of microscopic calculations^ is straightfor- 
wardly elucidated within the topological scenario of the 



5 



QCP. There are reports on the divergence of the effective 
mass in 2D hquid ^He, (see e.g. Refs. [sUslilli). Fur- 
thermore, following^, authors of Ref. [33 reported that 
at the density p > 9.00nm~^, the low-temperature hmit 
of the product Tx{T) quickly increases with increasing 
p. In addition, the ratio of the specific heat C{T) to 
T does not obey FL theory in this density region, since 
it increases with lowering T. These facts can be inter- 
preted as evidence for the presence of the FCi^ Unfortu- 
nately, so far the accuracy of extremely difficult mea- 
surements of the ratio C{T)/T at T < IK, is insuf- 
ficient to properly evaluate a low-temperature part of 
the entropy S and compare it with the respective FC 
entropy, extracted from data on x{T)- The divergence 
of the ratio C{T)/T, associated with the QCP, is ob- 
served in several heavy-fermion compounds.— i^Si^ Au- 
thors of the experimental article^- claim that data on 
the Sommerfeld- Wilson ratio Rsw = x{T)/C{T) in a 
doped compound YbRh2(Sio.95Geo.o5)2 point to an en- 
hancement of the Stoner factor that has to be infinite at 
the point of the ferromagnetic phase transition. However, 
evaluation of the Stoner factor from experimental data in 
heavy-fermion metals encounters difficulties, discussed 
in Ref. [271 . Furthermore, with a correct normalization 
experimental data^^ are explained without any enhance- 
ment of the Stoner factor^^ that rules out the relevance 
of the ferromagnetic phase transition to the QCP in this 
metal. A different H-T phase diagram is constructed for 



the heavy-fermion metal YbAgGe in Refs^TI. On this 
diagram, the FL phase is separated from a phase with 
magnetic ordering by a significant domain of NFL be- 
havior. Such a separation is easily explained within the 
topological scenario of the QCP, discussed in this article. 

In conclusion, in this article, two different scenarios 
of the quantum critical point (QCP), a low-temperature 
instability of the Landau state, related to the divergence 
of the density of states iV(0) ^ M* , are analyzed. We 
discuss shortcomings of the conventional scenario of the 
QCP, where the divergence of the effective mass M* is 
attributed to vanishing of the quasiparticle weight in the 
single-particle state. In a different, topological scenario, 
associated with the change of the topology of the Fermi 
surface at the QCP, the quasiparticle picture holds on 
both the sides of the QCP. This scenario is in agreement 
with microscopic calculations of the QCP in 2D electron 
gas and does not contradict relevant experimental data 
on 2D liquid "^He and heavy-fermion metals. 
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